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The practical focus of this work is the dynamical simulation of polarization transport processes in quantum 

spin microscopy and spectroscopy. The simulation framework is built-up progressively, beginning with state- 

^^ spaces (configuration manifolds) that are geometrically natural, introducing coordinates that are algebraically 

^~~^ natural; and finally specifying dynamical potentials that are physically natural; in each respect explicit 

^— ^ criteria are given for "naturality." The resulting framework encompasses Hamiltonian flow (both classical 

^ and quantum), quantum Lindbladian processes, and classical thermostatic processes. Lindbladian processes 

are shown to act gencrically to concentrate trajectories onto reduced-dimension state-spaces, such that the 

pullback of the symplcctic, metric, and complex structure of Hilbert space — the Kahler triple — is geometrically 

and informatically natural. The physical picture is that Lindbladian processes act to quench high-order 

correlations, such that the associated flow induces dimension-reduction. A concrete set of stochastic dynamical 

forms is provided, within both Ito and Stratonovich formalisms, that is Lindblad-complete and suited to 

numerical computation. These dynamical forms are shown to be the quantum counterpart of classical 

Langevin thermostats; they further encompass quantum measurement processes. Constructive validation and 

verification criteria are given for metric and symplectic fiows on classical, quantum, and hybrid state-spaces. 

"S Two illustrative classical cases are developed: Riemannian geodesic trajectories on a torus, and rigid-body 

^ rotational dynamics of water molecules. In particular, the projective naturality of a well-known quaternionic 

^ framework for simulating water molecule dynamics is proved as a special case of a general theorem that governs 

<sJ the construction of involutive subbundles upon canonically symplectic bundle manifolds. One quantum case is 

developed: dynamic nuclear-spin polarization (DNP) on tensor network state-spaces of varying rank. The 

,_^ resulting trajectory-oriented framework is well-suited to analyzing the practical performance and systems 

^ design challenges of classical, quantum, and hybrid systems like spin microscopes and spectrometers. 
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The work reported here has a practical focus: the end- 
to-end dynamical simulation of atomic-resolution quan- 
tum spin microscopy-iJ. Here "end-to-end" means an in- 
tegrated simulation of all the dynamical elements of a 
quantum spin microscope, from macroscopic elements like 
sample positioners, to mesoscopic elements like force mi- 
croscope cantilevers, to fully quantum elements, like the 
individual spins in supramolecular structures. 

Integrated simulations are necessary because present- 
day spin microscopes exhibit a level of complexity compa- 
rable to that of an earth-orbiting satellite. Like satellites, 
spin microscopes are multisystem devices that are chal- 
lenging to design, fabricate, operate, and repair. More- 
over, both satellites and spin microscopes seek to approach 
quantum and thermodynamic limits to performance as 
closely as feasible. 

To manage this complexity, it is routine engineering 
practice to simulate the performance of satellites (and 
comparably complex systems) end-to-end at every stage of 
design, fabrication, and operation. Our chief practical ob- 
jective in this article is to provide a comparable quantum 
systems engineering capability to spin microscopy. 

Among the most challenging aspects of systems engi- 
neering in spin microscopy is the quantum dynamics of the 
sample spins, which take place on a Hilbert state-space of 
enormous dimensionality. Even in the dynamically sim- 
pler context of magnetic resonance spectroscopy, existing 
quantum dynamical s imu lation methods typically scale 
poorly to large systemPEl or are "plagued by uncontrolled 
approximations"'^. Spin microscopy introduces new dy- 
namical elements — for example, ultra-strong magnetic 
field gradients — that further complicate the challenge of 
efficient quantum spin simulation. 

Quantum systems engineers regard classical molecular 
dynamics simulations with envy for their applicability to 
large, inhomogeneous dynamical systems. For example, 
it is routine practice to determine the diffusion constant 
D of liquid water by integrating the trajectories of many 
thousands of water molecules. The present article fo- 
cuses upon creating a simulation framework having a 
parallel capacity for quantum spin dynamics, specifically, 
having the capacity to determine spin polarization trans- 
port coefficients by numerically integrating the dynamical 
trajectories of hundreds-to-thousands of quantum spins. 
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A. Mathematical terminology 

To assist readers, we italicize each mathematical term 
as it is introduced, and we follow Terence Tao's exposi- 
tory practice of providing some explanatory context at 
the first appearance of each mathematical terrrP. For 
concision, though, we do not define those mathematical 
terms that appear in the set of canonical texts^ ^^ that 
we have listed in Sec. jIE] and we refrain almost entirely 
from introducing new mathematical terms.'^ For us, the 



^ first law of thermodynamics 

-» symplectic flows conserve energy 

?± Hamiltonian potential is Lie-invariant 

<- traverse ski-paths are constant-altitude 



second law of thermodynamics 

ensemble entropy is nondecreasing 
symplectic structure is Lie-invariant 
traverse skiers do not concentrate 




/\ 



metric flows 



dynamical 

conservation laws 

and thermodynamics 




dynamical 
concentration 
and pullback 



concentration & pullback 

metric flows concentrate trajectories 
metric structure is not Lie-invariant 
downhill skiers do concentrate 



reciprocity & causality* 

communication is reciprocal & causal 
Kahler-Lindblad-lto-Stratonovich pullback 
skiers mutually communicate & react 



'these principles are unique to quantum simulation 



FIG. 1. Geometric perspective of trajectory flow in classical and quantum simulations. Dynamical flows are specified in two 
ways: as How transverse to potential gradients (the symplectic flow of A-B) and as flow aligned with potential gradients 
(the metric flow of C-D); typically symplectic and metric flows are simultaneously present. By definition, a symplectic flow 
(A— B) specifies the trajectory tangent vector from a potential gradient via a symplectic structure. Symplectic flows are associated 
with Hamiltonian processes; they support conservation laws and respect thermodynamic principles. Metric flows (C-D) specify 
the trajectory tangent vector from a potential gradient via a metric structure. Metric flows are associated with thermostatic 
processes (both classical and quantum) and measurement processes (quantum only); metric flows concentrate trajectories and 
support pullback onto lower-dimension state-spaces. 



terms Hamiltonian and Lindbladian serve as both nouns 
and adjectives, and in particular the Lindbladian is a 
set of stochastic potentials (see Table IV I , originating 



in quantum information theory, that encompasses and 
generalizes the notion of a Hamiltonian. Our framework 
reduces simulation to the computation of trajectories as 
integral curves of Lindbladian potentials. 



B. Overview of dimension reduction 

As is common in the literature, in this article we take 
reduction to be a generic name for a constellation of 
mathematical tools that, given a dynamical system, "re- 
stricts the study of its dynamics to a system of smaller 
dimension. "iISl Particularly in the numerical simulation of 
quantum systems — for which the native Hilhert space has 
exponentially many dimensions — dimension-reduction is 
a practical necessity. 

Our approach to reduction focuses on thermodynamical 
processes; thus our framework is best suited to systems 
that are (or have been) in contact with a thermal reser- 
voir. In classical dynamics thermalizing processes are 
commonly called thermostats; in quantum dynamics they 
are commonly called Lindbladian processes. 

In this article we regard Lindbladian processes as be- 
ing fundamental, such that classical thermostatic pro- 
cesses emerge in a (well posed) quantum-to-classical limit. 
Hamiltonian potentials will govern symplectic dynamical 
flows, as usual in both classical and quantum physics, but 
even these Hamiltonian potentials will be subsumed into 
Lindbladian processes. 

Turning to a geometric point-of-view, state-space trajec- 
tories can be regarded as arising alternatively from metric 
structure or from symplectic structure; we will find that 
thermodynamical flows are metric, while Hamiltonian 



flows are symplectic. 

We therefore require of our simulations — both classical 
and quantum — that the state-spaces be endowed with 
both metric and symplectic structure. For thermodynamic 
processes in particular (and also for control processes) the 
dynamical potentials of our simulations will be stochastic^ 
and we will see that the stochastic nature of our potentials 
acts both to thermalize the dynamics and to enforce the 
standard quantum limits to measurement processes and 
backaction. 

We will find that Lindbladian reduction is well-suited 
to simulating large-scale quantum dynamical systems that 
lack symmetry (or in which the symmetry is present but 
unrecognized), but which are in contact with a thermal 
reservoir; quantum spin biomicroscopy is a paradigmatic 
example of this class of systems. 



C. Lindbladian versus Hamiltonian reduction 

Symmetries of Hamiltonian systems have long been 
known to enable dimensional reduction, and this has led 
to a constellation of closely related dimension reduction 
methods that are variously called Hamiltonian, symplec- 
tic, symmetric, geometric, or Lie reduction^^^^ . These 
various names reflect, not varying methods, but varying 
emphasis: "Hamiltonian" as a modifler emphasizes the 
role of dynamical potentials, "symplectic" emphasizes that 
the dynamical flow is a symplectomorphism, "symmetric" 
emphasizes the central role of symmetry in the dimension 
reduction, "geometric" emphasizes considerations relating 
to differential geometry, and "Lie" emphasizes the central 
role that Lie groups play in these methods. 

Our simulation framework comprises an alternative 
approach to dimension reduction that can variously be 
called Lindbladian, metric, thermostatic, or synoptic. Here 



"Lindbladian" emphasizes the key role of quantum stochas- 
tic processes, "metric" emphasizes the role of Riemannian 
metric flow, "thermostatic" emphasizes reduction arising 
from dynamical interaction with a thermal bath, and "syn- 
optic" emphasizes that Lindbladian processes associated 
with thermalization can always be described as equiva- 
lent continuous measurement processes that synoptically 
observe the dynamics^"^. 

Our preferred phrase will be "Lindbladian reduction" 
in recognition of the key physical idea is that dynamical 
flows associated with noise, thermalization, or measure- 
ment (Lindbladian processes accommodate all three) act 
generically to compress trajectories — both classical and 
quantum — onto state-spaces of lower dimension. 



D. New results in this article 

We now list those of our results that are new (to the 
best of our knowledge), together with the sections in 
which they are derived. 

Lindbladian reduction The expressions given in Ap- 
pendix [A] and in Table |IV| express general Lindbladian 
processes in terms of stochastic flows that are specified 
in terms of potentials known as Berezin symbols. 

Compatible structures Theorems 1-3 of SecHlTland 
Tables [lf[Tl| link metric structure compatibly with sym- 
plectic structure. 

Block potentials The Block equations provide a 
much-used paradigmatic description of spin thermaliza- 
tion; in Appendix |A] we derive a complete set of Bloch- 
equivalent Lindbladian potentials. 

Factored metric The metric structure specified in 

Table 

is ena 



inrj-L) is given in a factored representation that 
Dling for large-scale simulation codes. 



Validation and verification criteria Large-scale 
simulation codes are themselves engineered objects; the 
validation and verification criteria of Table ^E) and 
Table ^JJF) are tools for creating these codes. 



E. Dovetailed elements of naturality 

We have organized our framework around dovetailed 
notions of mathematical naturality. For convenience we 
have identifoed a set of twelve canonical textbooka^Hi^ 
that, in aggregate, convey all of the elements of naturality 
that we require. We review them in a recommended 
reading order, which we also adopt as a descending order 
of precedence for mathematical notation. Thus the ideas 
of mathematical naturality that are most important to 
our simulation framework are systematically given the 
highest notational precedence. 

Geometric naturality We take geometric naturality 




FIG. 2. Geodesic trajectories on a torus. The trajectories are 
numerically integrated upon an algebraic state-space that is 
constructed via Theorem 1 (Compatible Metric Structure) and 
Theorem 2 (Compatible Symplectic Structure); see Tables [l[|ll| 
and Sec. IIIB (A-B) convergent geodesic trajectories mainly 
traverse regions of positive Riemannian curvature, while di- 
vergent geodesic trajectories (C-D) mainly traverse regions 
of negative curvature (see Lee [7j Ch. 10 for details). The 
starting orientation of trajectory (C) differs from those of (D) 
by ~ 1.0x10"^'' radians; both trajectories circle the central 
hole six times before returning near to their starting point; 
the (C) trajectories thread the central hole and the(D) tra- 
jectories don't. Algebraic manifolds are well-suited to the 
high-accuracy, numerically efficient trajectory integration that 
is required to discriminate these fine differences. 



to have the highest mathematical precedence, so that 
the elements of our simulation framework are specified in 
geometric terms whenever possible. The key elements of 
geometric naturality are taken to be the pullback of dy- 
namical potentials to the simulation state-space, the spec- 
ification of compatible metric and symplectic structure, 
and the pushforward of trajectories to the native state- 
space. We adopt the notation and terminology of Lee's 
Sm,ooth Manifold^ and Riemannian Manifold^, with 
minimal extensions to complex manifolds as described in 
MartirP. With one necessary exception (see the heading 
"algebraic naturality" below) the elements of our simula- 
tion framework are compatible with geometrically natural 
constructions; this allows the development of simulation 
software to be largely automated. 

Dynamical naturality The key elements of dynami- 
cal naturality we take to be a symplectic flow governed 
by Hamiltonian potentials, and metric flow that describes 
thermalizing dynamics and measurement processes. We 
adopt the notation and terminology of ArnolcP and of 
Abraham and MarsderpSl^ with extensions for classical 
thermostatic dynamics from Frenkel and Smit"'^^. Mag- 
netic resonance dynamics, and in particular Bloch relax- 



ation, are specified in the notation of SlichteplS. Dynami- 
cal naturality is integrated with geometric naturahty by 
the geometrically natural pullback of dynamical poten- 
tials, metric structure, and symplectic structure. 

Projective naturality Combining the preceding no- 
tions of geometric and dynamical naturality leads us to a 
notion of projective naturality that is concretely realized 
as the validation and verification criteria of Table lltE) 
and Table [ll|^F) , which we construct by methods that are 
geometrically natural. The key role in systems engineer- 
ing of validation and verific ation criteria are described 
in a classic article of BoehnP^IM projective naturality 
amounts to the expression of Boehm's criteria in a form 
that is geometrically and dynamically natural. 

Informatic naturality In our framework all pro- 
cesses of thermalization, measurement, and control are 
realized as stochastic Lindbladian flows, which we develop 
by merging the quantum informatic approach of Nielsen 
and Chuang (il4t chs. 2 and 8) with the stochastic ap- 
proach of Kloeden and PlaterfiS Nielsen and Chuang 
specify Lindbladian dynamics in an algebraic form that 
our simulation framework translatse into one-forms speci- 
fied by Berezin symbols (see Appendix [A| for recipes; see 
(P51 sec. 3 and figs. 8-9) for physical systems that exhibit 
Lindbladian dynamics). The resulting general Lindbla- 
dian is natural in geometric, dynamical, and projective 
senses described above; this too helps in automating the 
development of simulation software. 



Algebraic naturality For reasons of simulation ef- 
ficience, we design the symplectic and metric structure of 
our simulations to have an algebraic structure that we as- 
sociate to a choice of algebraic coordinates (see Sec. Ill C ). 



The specification of algebraic coordinates is the sole as- 
pect of our framework that is not natural the geometric, 
dynamical, projective, and informatic senses described 
above. Once coordinates have been specified, the ensu- 
ing trajectory integration is carried through by standard 
numerical computation methods that are adapted largely 
from Golub and Van LoarP^ and from Greenbaurrfl^ 

Further aspects of naturality Our framework fo- 
cuses on the simulation of noisy, asymmetric systems; we 
therefore have little occasion to refer to the literature on 
symmetric reduction methoda^^lEI]^ despite the undoubted 
utility of symmetry reduction in particular cases. 

Similarly, our quantum simulation state-spaces in gen- 
eral are not Hilbert spaces, but rather are lower-dimension 
nonlinear tensor network manifolds. Thus we will have 
little occasion to apply such familiar tools of quantum me- 
chanics as density matrices or projective measurements, 
because these are not defined on tensor network manifolds. 

Our own recent analysis^Sl provides derivations of many 
of the standard results of quantum mechanics — including 
uncertainty principles, standard quantum limits to noise 
and measurement accuracy, and projective measurements — 
upon this non-Hilbert quantum state-space. This article 




FIG. 3. Rigid body rotation of a water molecule, accompa- 
nied by center-of-mass motion in a harmonic potential. As 
with the toric trajectories of Fig. [2] this motion is formally 
a geodesic trajectory upon an algebraic state-space that is 
constructed via Theorems 1-3 of Tables UHlll The rotational 
sector of the state-space so constructed is essentially ide ntical 
to the quaternionic state-space of Miller et alf^, see Sec. HIE 
for details. The molecule's angular momentum is initialized 
to be nearly aligned with the H2O molecule's middle moment 
of rotational inertia; in consequence the rotational motion 
consists of intervals of stable precession alternating with in- 
tervals of rotational tumbling. Astronaut Michael Foale has 
given a diverting account of the practical consequences of this 
long-known tumbling instability^. 



is recommended to readers who wish to work through the 
details of how results typically derived on linear Hilbert 
spaces can be alternatively derived via a Lindbladian 
formalism on tensor network manifolds. 

Educational objectives It plausible that during the 
last two centuries, more people have learned the ba- 
sic principles of differential geometry from Nathaniel 
Bowditch's nautical textbook The American Practical 
Navigator (commonly known as ^^Bowditch") than from 
any other single book or article. It is less widely appre- 
ciated, however, that editions of Bowditch as early as 
1807 are startlingly rigorous in their mathematical de- 
scription of what became known as Riemannian geometry 
d^SJ p. 91), decades in advance of Gauss' pioneering work. 
What accounts for the enduring success of Bowditch 
over the past two centuries? There is a storjPSl that 

. . . Bowditch vowed while writing (the first 
American) edition to "put down in the book 
nothing I can't teach the crew," and it is said 
that every member of his crew including the 
cook could take a lunar observation and plot 
the ship's position. 

Thus Bowditch provides a vivid historical example of Mac 
Lane's point that "Mechanics developed by the treat- 
ment of many specific problems" (P71 p. 295) and the 
observation of Abraham and Marsden that "In mechanics 



almost every branch of mathematics gets used" (jTUl p. 15). 
The three worked examples in this article were crafted, 
therefore, with conscious attention to the principles of 
Bowditch, Mac Lane, Abraham, and Marsden. 

Scientific objectives As we gain a broader appreci- 
ation that the dynamical elements of quantum mechanics 
survive puUback onto nonlinear state-spaces, the natural 
question arises "Is the state-space of quantum mechanics 
a linear Hilbert space, or does it have some more general — 
even dynamical — geometry?" Formally speaking, in this 
article we view the state-space of nature as a Hilbert space, 
to which our simulation framework provides a useful ap- 
proximation. There is a substantial body of literature, 
however, that takes the opposite point-of-view, which we 
review in posing Challenge |4] of Sec. |VD[ 



II. THE ASCENT: GEOMETRIC DYNAMICS 

We start with the idea that Hamiltonian simulation 
frameworks — whether classical or quantum — can be con- 
structed from objects that live naturally on non-flat state- 
spaces. Our starting inventory is curves and real- valued 
functions. From curves we construct tangent vectors and 
from functions we construct gradient one-forms. Then 
by taking products we construct symplectic forms and 
metric forms; these bilinear forms define our state-space 
geometry. 

The classical Hamiltonian framework of Arnol'cP' illumi- 
nates the dynamical roles of symplectic and metric forms 
as follows. From the gradient one-form of a Hamiltonian 
potential a vector field is constructed that is called the 
Hamiltonian flow. Then simulation trajectories are simply 
the integral curves of the flow (Fig. nK-B) . 



A. Roles of symplectic flow 

The symplectic form respects the first law of thermody- 
namics (energy conservation) by answering the question, 
"How fast and in what direction do we move?" as follows: 
"Given the gradient one-form of a Hamiltonian potential, 
here is a vector field for which that gradient vanishes; 
flow tangent to this vector field." Figure [l]^ provides the 
analogy of skiers traversing a mountainside, with altitude 
in the role of a Hamiltonian potential. 

The second law (nondecreasing entropy) is respected by 
requiring that the symplectic form be closed, so that its 
Lie derivative vanishes, which is the geometric statement 
of Liouville 's theorem. Physically this means that under 
symplectic flow, state-space density does not concentrate. 
Figure [Tj3 provides the analogy of a crowd of skiers whose 
density is unchanging. 

Mathematically speaking, symplectic forms specify a 
natural symplectic isomorphism between one-forms (in 
this case, the gradient of the Hamiltonian potential) and 
vector flelds (the Hamiltonian flow). 
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FIG. 4. We suppose that on a native state-space manifold A'^, 
the dynamical potentials and the state-space metric g^ have 
been specified, and we further suppose that the simulation map 
F: M ^ N has been specified. Considering first the native 
state-space, the following are geometrically natural construc- 
tions for a complete set of tangent and cotangent subbundles: 
T*N\j, = KerF*, TN\j, = ImK, TN\j, = 0;^ {T*MM = 

{TNM^, and r*iV|x = <l>g„ {TMM = {T*NM^. In con- 
trast, on the simulation state-space only two of the four sub- 
bundles have natural constructions, namely T*M\ji = ImF* 
and TM\yi = Ker_F*. The remaining unspecified elements are 
labeled "?"; these must be designed, validated, and verified 
in service of engineering objectives, via the constructions and 
theorems of Section Hill 



B. Roles of metric flow 



Now we turn our attention to flows that are metric 
rather than symplectic. Classical protein simulations 
often simulate flows induced by interaction with a low- 
temperature thermal bath, to find (by simulation) protein 
folds that are low-energy and hence stable. At physiologi- 
cal temperatures covalent bonds are effectively rigid links; 
the conformational state-space then is parameterized by 
bond angles and has lower dimension than the Cartesian 
space in which it is immersed. The pulled-back Cartesian 
metric endows the bond-angle state-space with a metric 
structure, and by setting the kinetic terms in the Hamil- 
tonian to zero, the state-space is further endowed with a 
chemical energy function. 

These two new structures — the metric form and the 
energy function — provide a new answer to the question 
"How fast and in what direction do we move?" , namely, 
"Here is a vector field that maximizes the energy function's 
rate-of-change; flow tangent to this vector field" (Fig. [l]C). 
Like symplectic forms, metric forms establish a metric 
isomorphism between one-forms (the energy gradient) 
and vector fields (the flow of steepest descent) . 

We thus arrive at a classical simulation framework 
whose elements are geometrically natural, in which tra- 
jectories are governed by symplectic and metric flows 
specified in terms of potentials and gradients, and we see 
how to gain efficiency by adapting state-space geometry 
to trajectories. 



C. Quantum symplectic flow 

To extend these principles to quantum systems we first 
review the symplectic elements of quantum dynamics. 
These are well understoocP^H^ni^ and we need only specify 
them in a form that is suited to our geometric simulation 
framework. 

Hamiltonian trajectories — both classical and quantum — 
are the integral curves of the Hamiltonian flow v (a 
vector field) that satisfies (/'tj(w) = ViJ, where H is the 
Hamiltonian potential, V7J = dH is a gradient one- 
form, w is a symplectic form, and (in various notations) 
4>u>{v) = w(u, ■) = iv^ = V juj is the one-form that is in- 
duced jointly by uj and v. 

Some authors conventionally define a vector-valued 
symplectic gradient to be V(^ = ^J^'^V, such that Hamil- 
ton's equation is w = V^^iJ; we prefer forms because forms 
puUback naturally and vectors don't. 

To link-up symplectic dynamics with quantum physics, 
the Schrodinger equation 4>g{v) — ~ J*Vi?/(2fi) must be 
satisfied, where h is Planck's constant, G is the Hilhert 
space metric, 4>g{v) = G{v, ■) is the metric isomorphism, 
and J* is a complex structure on one-forms (that is, an 
automorphism satisfying {J*Y^ = — </*; thus for us the 
sym,plectic gradient J*V maps functions — >■ one- forms). 

Consistency requires that the Hamilton and Schrodinger 
equations predict the same trajectory. By direct substitu- 
tion this requires that the bijective maps {4>G,4'unJ*} be 
a compatible triple under composition " o " , i. e. (f>G— — 
J*o0„/(2fi), (1)^ = 2fiJ*o^G, and J* = (p^o(p-^^ /{2h). 
This triple induces a complex structure on vectors 
J = (f)Q o J* o (j)Q that is a symplectic and metric isome- 
try, G{J-, J-) = G{-, •) and uj{J-, J-) =w(-, •) (as is usual 
in quantum mechanics). 

Geometers often further require that J be integrable, 
which we ensure by inducing J via holomorphic coordi- 
nates; the Newlander-Nirenberg theorenPil then guaran- 
tees that J is integrable. The practical advantage is that 
the natural action of J (J*) on the vector (covector) bases 
that are induced by holomorphic coordinates is simply 
multiplication by factors of ±i. 

The physical lesson is that Schrodinger dynamics is 
symplectic — and thus thermodynamical — if and only if 
any two of the symplectic, metric, and complex structures 
{uj, G, J} compatibly determine the third. A state-space 
that is endowed with an integrable compatible triple is, 
by definition, a Kdhler manifold^^; specifying quantum 
dynamics as a symplectic flow on a Kahler manifold is 
called geometric quantum mechanicP-^-, this concludes 
our summary of it. 



D. Quantum metric flow 

Symplectic flow alone is insufficient for quantum 
simulation — ^just as it is classically insufficient — and so 
we turn our attention to quantum metric flows. Metric 
flows appear naturally in processes involving noise and 



measurement. These processes are called quantum opera- 
tions, and their well-known ma them atical expression is 
the Lindbladian master equatiorP^^. 

A central theme of our geometric simulation frameworks 
is that when the Lindbladian equation is unravellea^"^ in 
terms of Hamiltonian potentials and their gradients, the 
resulting flow is concentrative. 

This leads us to envision a quantum simulation frame- 
work (Fig.[lp-D) that has the same geometric structure as 
classical simulation, in which symplectic flow (Schrodinger 
rather than Hamiltonian) unites with metric flow (Lind- 
bladian rather than thermodynamic) and is pulled-back 
onto a lower-dimension state-space (Kahlerian rather than 
bond- angle) . 

Quantum master equati ons ca n be unravelled in var- 
ious gauge-equivalent way^^^'^, of which the synoptic 
gaug^^ is well-suited to implementing simulations. This 
gauge exploits a reciprocal informatic symmetry^ "We ob- 
serve a system only by allowing the system to observe 
us." Physically speaking, the synoptic gauge simulates all 
noise processes as equivalent measurement processes; in ef- 
fect this describes quantum operations as communication 
channels. 

Our skiing metaphor therefore depicts each Lindblad- 
Ito skier as descending rather than traversing a mountain- 
side, while observing and responding to the motions of 
other skiers . . . who in turn observe and respond to other 
skiers . . . etc. We will find that Lindblad-Ito one-forms 
have a stochastic component that reflects this mutual 
observation-and-response. 

During the descent we will give detailed account of 
quantum Lindbladian processes, but flrst we survey the 
geometric summit at which we have arrived. 



III. THE SUMMIT: VALID SIMULATION DESIGN 

Our ascent has been abstract and rapid — now we survey 
the summit with a view toward planning our descent to 
practical applications. Here we are motivated by Ed Vi- 
estur's mountaineering maxim p5l p. 168) that "Getting 
to the top is optional. Getting down is mandatory." 

To assist us with planning our descents to practical 
applications, we are going to prove three powerful but 
abstract theorems that relate simulation verification and 
validation (Tables [iHh| . We motivate these proofs in the 
light of two practical examples: geodesic flow on a torus 
(shown in Fig. ^ and rigid body rotation of a water 
molecule (shown in Fig. p| 



A. Reconciling with Grotfiendieck 

In discussing our example, we bear in mind a famous 
passage of Grothendieck that describes his ideal of math- 
ematical naturality: 

The unknown thing to be known appeared 
to me as some stretch of earth or hard marl. 



resisting penetration. . . . the sea advances 
insensibly in silence, nothing seems to happen, 
nothing moves, the water is so far off you 
hardly hear it. ... yet it finally surrounds the 
resistant substance. 

Regrettably, our experience as educators has been that 
very many students in engineering and science find a 
Grothendieck-style mode of thinking to be so abstract 
as to be dauntingly alien; these students learn more eas- 
ily in the context suggested by Mac Lane's aforemen- 
tioned maxim "Mechanics developed by the treatment 
of many specific problems" (|27l p. 295) and Arnold's 
aforementioned "Hamiltonian mechanics is geometry in 
phase space; phase space has the structure of a symplectic 
manifold" ^ p. 161). 

To reconcile these various styles of naturality, with a 
view toward making natural simulation frameworks easier 
to grasp, we are going to motivate our natural-style theo- 
rems by Arnold-style and Mac Lane-style consideration 
of practical examples. In particular, for us the sea of sim- 
ulation naturality will not "advance insensibly in silence" 
until it "surrounds the resistant substance." Rather, we 
will take the point of view that naturality is itself a utili- 
tarian construct of engineering and science, and we will 
motivate the development of our naturality theorems "on 
the fly" with respect to concrete aspects of our sample 
problems. 



B. Introducing coordinates 

We begin with the intuition that the torus is a two- 
dimensional surface embedded in a three-dimensional Eu- 
clidean state-space N. We take Cartesian coordinate func- 
tions N to be {x^, x'^, x^}, and we take coordinate func- 
tions on the simulation state-space M to be {(j°, q^,q^,q'^}, 
which are taken to be related to {x^, x'^, x^} by 



rxM 




x^ 


= 


x^ 





{ri + ^2 cos 9) cos (f) 






{ti + T2 cos d) sin 






r-i sin Q 


= 


arctan(go, Si) 




e = 


arctan(g2,g3) 



(1) 



Already we have departed substantially from traditional 
descriptions of toric geodesies, which are based upon 
the above Clairaut coordinates {9,(j)} (see p. 219 of [55)) , 
in that we have introduced four coordinate functions 
{q^ T q^ , q^ , q'^} to simulate motion on a two-dimensional 
manifold. Thus the map F : M —^ N \s now specified, 
but it has not been specified in a geometrically natural 
manner. Rather, by designating an overcomplete set of 
coordinates {<z°,q^,'Z^, 9'^} we will obtain an algebraically 
efficient simulation. We will see that choosing coordinates 
for F is the only step of our simulation framework that is 
not determined entirely by considerations of naturality. 
We further associate to the Euclidean coordinates 
{a;^, a;^,a;'^} the momenta {p^,p^,p^} and the Hamilto- 
nian function H = ((p^)^ + {p^Y + {p^Y) /(2m). Requir- 
ing that the Hamiltonian be a scalar associates to our 



TABLE I. Providing a compatible metric structure to simu- 
lation state-space. See Sec. |III| for physical motivation and 
notational details. 

A A smooth map F: i\/ — )• Af from a simulation state-space M to 
a native state-space N induces a natural pullback of the thermo- 
dynamic potentials that govern metric flow: 

't'afS--^^) — ^Gn (metric flow is determined by a 

Qj^j ^ p o Qj^ thermodynamic potential Gjv on N 



Kri^M) = VG, 



whose pullback to M is natural) 



B By definition, a metric flow on M is called F-compatible if the 
pushforward of the simulated metric flow Xm is a projection of 
the native state-space metric flow Xn- 



F,{Xm) 

TTe^ O TTe^ 



T^F, 



(a metric g^i that induces an 
_F-compatible flow is said to be 
an F-compatible metric) 



C Theorem 1 (Criteria for a Compatible Simulation Metric) 

Given the natural subbundles r*i\/|3j* = luiF* and Ti\/|^ = 
KerF,, choose T*Af|K* and TJ\/|;r to be any complemen- 
tary subbundles, that is, TM = TM\ji®TM\j^ and T*M = 
T*M|k*®T*M|n*. Then gu is F-compatible iff: 

D A diagram-chasing proof of Theorem 1 (maps labeled with a "~" 
are isomorphisms): F, 



§« 


TMU 


t^ 


X 


^-: a. 






TO u^ 


K 


2 


=1 <{> 


^ 13 




t^ 


-^ 


,' 


'"-a 


T*A 


/|.r 



TM\^ ; 


TWU 


. 


t 






I 1 


^ 


; 


N 


1 


^»., 


T'M\^. ' 


T*WU 



TTVI, 






u 
■^ I/) 



A set of projectors for verifying that simulations of metric flow 
are F-compatible: 



TTF' 



F. o ,^-1 O F* O ^, 



,Pg^ O F-, O (j,- 

F* o <t)„ 

't>'al O F' O 4>„ 



O F* 

° 4>a 
^9m 

o F, 



(the 7r- each pass 

the verifying test 

of idempotency: 



F By definition, an F-compatible metric gM is called F-natural if 
a coordinate basis exists on M (also called a holonomic basis) for 
which the duals of the natural subbundles T*A/|j;» and TAI\t^- 
are orthogonal complements. We may then choose TMIj; = 
{T'MM' = (TMM^ and T*M\y,. = {TMM* = (T'MM^ 
(for which the two dashed-diagonal lines in (D) both vanish). 

native state-space a natural metric g^ ■ TN -^ T*N that 
evidently is given by 



5JV = 2m [dx^ ® dx^ + dx^ ® dx"^ -f dx^ 



• dx-" 



(2) 



Here we are including the mass m because when extend 
this construction to the multi-atom water molecule of 
Fig. |3l the individual atoms will have varying masses. 

From now on, and in all our examples, we will regard 



the native metric gisi and the simulation map F to be 
specified ab initio. 



TABLE II. Providing a compatible symplectic structure to 
simulation state-space. See Sec. |III| for physical motivation 
and notational details. 



C. Constructing natural simulation metrics 

It is evident from the many "?" markings in Fig.ldlthat 
the preceding construction is insufficient for specifying 
a complete set of natural structures on our simulation 
state-space. The missing subbundle structures must be 
designed, validated, and verified in service of engineering 
objectives, and specifying that complete set of natural 
structures is the objective of (Tables IMl]) . 

We require that our simulations be natural for both 
metric and symplectic dynamics, and we consider the 
(simpler) metric case first, which is given in Table ^ The 
notation of Table |T] follows without essential change the 
notation of Lee (jH see pp. 130 and 317-19). We write 
(lA) to reference the first entry in this table; we adopt 
this concise notation for all this article's (many) table 
entries, from (1A)-(4I). 

Our key engineering objective is that our simulated 
trajectories be projectively natural in the metric sense of 
(IB) and in the symplectic sense of (2B). 

Given a native metric g^ and a simulation map F, (IF) 
provides criteria for endowing the simulation state-space 
with a metric gM that is F -natural. Constructing an 
_F-natural g^ is the key step in specifying projectively 
natural simulations, since all the missing elements labeled 
"?" in Fig. |4] can be filled-in once gM has been constructed. 

The following is a concrete construction of an F-natural 
simulation metric for the torus trajectories of Fig. [2] The 
sole inputs to the construction are the simulation map 
F of (fTl) and the native state-space metric gN of (l2]). 
The construction is can be carried through by automated 
programgSZl. 

In accord with (IB) the pulled-back native metric F* o 
gN ° F^, specifies two natural subbundles on M, namely 
T*M\ji* and TM\j^. Although the subbundles themselves 
are natural objects, there is in general no unique natural 
choice of basis vectors for them. By a singular value 
decomposition (SVD) of the matrix representation of F* o 
gN o F^ in the chosen coordinate basis, we may choose 
the subbundles to be spanned by the SVD basis vectors: 



T*M\ji, = Span {q^dq° - q"dq\ 
q^dq^ -- q^dq^ ) 

TM\yi = Span {q^d^o + q^d^i, 
q^dg2 + q^dqz ) 



(3) 



(4) 



Now comes the sole step in the construction that is not 
geometrically natural: the complementary subbundles 
TM\ji* and T*M\j4* are specified as coordinate duals: 

TM\ji = (T*M\ji.)* and T*M\y.i. = (TA/U)*. The sub- 
bundles so constructed then satisfy 

E e TAf Ik, and (5) 
, E e TM\oi (6) 



: for e G T*M\ji* 
e{E) = for e e T*Af |j^. 



e{E) 



A A bundle map F: M —>■ N from a simulation cotangent bundle 
M to a native cotangent bundle N induces a natural pullback of 
the Hamiltonian potentials of canonical symplectic flow: 






(symplectic flow is determined by 

a Hamiltonian potential if^y on N 

whose pullback to M is natural) 



F.(XM)=^^.(Xiv) 

TTr O TTfi 



B By definition, a symplectic flow on the bundle manifold M is said 
to be F-compatible if the pushforward of the simulated symplec- 
tic flow Xj^i is a projection of the native state-space flow Xf^: 

(a symplectic structure ujj^j that 
induces _F-compatible flow is said 
' E, ^ " E, ~ " F, to be an F-compatible structure) 

C Theorem 2 (Natural Metrics Induce a Foliation Isometry) 

Let TA/|;r = (Ker F,)-^ be the subbundle naturally associated 
to an _F-natural metric; then TJ\/|;r is an involution that induces 
a foliation of M for which F is an isometry on the leaves. 

D Theorem 3 (Construction of a Compatible Bundle Map) 

An F-natural metric gu induces a bundle map F such that the 
canonical symplectic structure lj^j on M is F-compatible: 



F|s: M - 



N : 



Fla-: T*M ^T*N 






(Facts on the leaves of the 
F-natural foliation of A/ 
as a symplectomorphism) 



E A diagram-chasing proof of Theorem 2 (maps labeled with a 
"~" are isomorphisms, and the dotted lines vanish in F-adapted 
normal coordinates): 

F.|:r^^ F,U^^ 







T*JV|j,- 



FI.3.-.3. 



F A set of projectors for verifying that simulations of symplectic 
flow are F-compatible: 



TTp = Ft o 6^ o F* o 6 

TTp. = F* o 6^ o F» o d)^^ 
71,-1 = cbT.'^ o F* o <A„ o Ft 



(the 7r each pass 

the verifying test 

of idempotency: 



Note that the above subbundle relations are geometrically 
natural, and thus hold in all coordinate systems, even 
though they refer to subbundles that were constructed by 
SVD methods in a particular coordinate system associated 
with (IllQ. Because these relations hold in all coordinate 
systems, the remainder of our analysis can be carried 
through entirely via natural constructions. 

Following the prescription of (IC) we then construct 
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an F-natural metric qm having the property that the 
above subbundles are simultaneously coordinate duals 
and orthogonal complements, as specified in (IF). Such 
metrics are said to be F -natural because, if we construct 
(suitably) Riemannian normal coordinates on M and N, 
then it is straightforward to verify that the four bundle 
maps {0gjy , 0gj,j^ , F^,, F*} can be made locally diagonal by 
a suitable adaptation of these coordinates. We will not 
give the details of the normal coordinate construction, 
which is straightforward but rather lengthy; the required 
computational idioms are thoroughly covered in Lee's two 
books (see Thms. 7.8-13, 11.24, 19.10, and A.33 of ref. IS 
or alternatively see Props. 5.7 to 5.11 of ref. [7|)- 

In the adapted normal coordinates so constructed. The- 
orems 2 and 3 of Table |TT] are trivially true, and because 
the theorems are covariantly stated, they hold in all coor- 
dinate systems; similarly the verification operators tt^ are 
(2E) projectors relative to any basis. Thus the main utility 
of normal coordinates is theorem-proving; in integrating 
dynamics trajectories we need not compute them. 



D. Dimension reduction and augmentation 

In essence what the theorems and constructions of 
Tables |l]-[I] provide is a general method, not only for di- 
mension reduction, but also for dimension augmentation. 
The theorems and constructions are particularly useful in 
cases in which dimension augmentation allows dynamical 
equations to be cast into a computationally advantageous 
algebraic form. 

As we will discuss later on, in the context of quan- 
tum simulation, Hilbert space can be be viewed as a 
dimensional augmentation of the (unknown) underlying 
quantum state-space of nature; thus dimension augmen- 
tation methods are of considerable fundamental interest 
in themselves, in addition to their practical utility in 
simulations. 



E. Rigid body rotation 

Now let us construct a dynamical simulation of rigid 
body motion by the methods of the preceding section. 
Our example will be the tumbling water molecule of Fig. [3) 
and the symplectic dynamical equations that we derive 
will be of a quaternionic variety has found much use in 
astronomy2SJ and in molecular simulation^. Our main 
goal here is to show that once quaternionic coordinates 
have been specified, then the derivation of the equations 
can be carried through automatically, via Theorems 1-3 
of Tables III ITTl and their associated natural constructions. 

Let Euclidean coordinates on the native state-space be a 
set of triplets {{xj , x^ , xf)}]i € {1, . . . , riatom) associated 
to canonical momenta {{pj,pf,p^)};i £ (1, . . . , natom)- 



Let the Hamiltonian potential be 



H = 



"atom 

E 

1=1 



{ipir+iph' + ipff)/^^^) (7) 



where mi is the mass of the i-th atom of the rigid body. 
We require that this potential be a geometric scalar, and 
from this requirement we read-off the dynamical metric 
gN of the native atomic state-space to be: 



9n 



E 2^i2m,){dxi<i^dxl) 
i=i j=i 



(8) 



We specify the coordinate map F by 



[i?(9°<z\'z','Z^): 



r(^').i 




\Q'^ 


i^'h 


— 


Q' 


[ix%\ 




W\ 



\Xf))i 



(9) 



where by inspection {Qi,Q2iQ^) are center-of-mass coor- 
dinates and i?(q°, g^ g^, q^) is a 3 x 3 rotation matrix that 
is a function of quaternionic coordinates {g", g\ g^, q^}. 

A quick way to derive the quaternionic functional form 
of i?((7°, q^, q^, q^) is via Pauli spin matrices {cr", ct^, a"^, cr"^} 
in which case the rotational transform R{q) : x ^^ x' \s 
given by 



(i'-^) = i?^(q)(x.rT)i?(q) 



where 



R{q) 



q"a" 



iq^a^ 



iq^a^ 



iq^a^ 



;(g0)2 + (gl)2 + (92)2 + (q3)2). 



(10) 



(11) 



We recover the standard quaternionic representation for 



R acting on the unit vector x by taking traces of ( 10 1 
with (11) substituted. 



Then just as we calculated the two natural toroidal 
subbundles T*M\ji* and TM\j^ of ^^, we calculate the 
two natural rigid body subbundles to be 

T*A/|3j. = Span( dQ^,dQ^, dQ^, (12) 

- q^ dq° + q" dq^ - q^ dq^ + q^ dq^ , 

- q^ dq° + q^ dq^ + q° dq^ - q^ dq^ , 

- q^ dq^ - q^ dq^ + q^ dq^ + g" dq^) 

TMU = Span {q^^d^o + q^d^i + q^dq2 + q^dq^) (13) 

The simple algebraic form of these quaternionic subbun- 
dles is, from a natural point of view, the sole (and suffi- 
cient) justification for embracing quaternionic coordinates. 
The remainder of the rigid body dynamical calculations 
shown in Fig. [3] proceeds according to the theorems and 
constructions of Tables [iHlTl precisely as for the previous 
case of geodesic motion on a torus. We will not give 
details, because the successive steps each are geometrically 
natural and the expressions obtained are equivalent to (for 
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example) the NO_SQUISH formalisir!^ of the molecular 
dynam.ics literature. 

The chief practical advantage of deriving the rigid- 
body dynamical equations by this route is conceptual 
simplicity: Theorems 1-3 guarantee that the simulation 
dynamics are symplectic, the proofs of the theorems are 
diagrammatic, the associated constructions are generically 
natural and thus can be automated, and the concrete 
validation and verification criteria of (IF) can be used to 
check the resulting dynamical simulation codes. 

It is no use pretending, though, that natural simulation 
methods inevitably produce dynamical simulation codes 
that are particularly short or simple. Even for our simple 
example of a water molecule, the projector tt^^^ of (IF) is 
an 18 X 18 matrix whose associated projective verification 
TTi^iv ■ I'l.cjjv = ""i^N requires the separate cancellation of 
18^ = 324 matrix entries. 

The simulation of rigid body motion in quaternionic 
coordinates thus marks a level of software complexity at 
which the availability of concrete methods for validation 
and verification starts to become a practical necessity. 



B. Lindblad-lto one- forms 



In (3C) we give the Lindblad-ItdSMEol one-forms asso- 
ciated with measurement of a general operator q, which 
are seen to have both a drift component and a diffusion 
component. These one-forms are given in the synopticP^l 
Lindbladian gauge that concentrates trajectories. We see 
that only the metric-flow gradient V enters in the one- 
forms — not the symplectic-flow gradient J*V — in accord 
with the downhill skiing metaphor of Fig.[lp-D. 

The Lindblad-lto dynamical forms of (3C) supply us 
also with a continuous measurement of (/"(t); this provides 
causal Hamiltonian control via (3D). 

In the appended material the dynamical forms (3A-D) 
are elaborated into recipes for simulating arbitrary Lind- 
bladian processes (a detailed derivation is given). The 
resulting dynamical forms are simple variations on (3A-D), 
with Block relaxation as a special case. Thus (3A-D) com- 
prise a complete set of quantum postulates — equivalent 
to Postulates 1-4 of Nielsen and ChuangSl^ for example — 
given in terms of geometric forms and dynamical flow. 



IV. THE DESCENT: QUANTUM SIMULATION 

Now we have gathered all of the physical ideas and 
mathematical tools needed to construct a quantum simu- 
lation framework in three explicit stages: flrst we spec- 
ify orthodox quantum dynamics in terms of forms on 
Hilbert state-spaces; then we pull-back the forms onto low- 
dimension Kahlerian state-spaces; and finally we specify 
efficient algorithms for integrating the flow of trajectories. 
In detail, this construction is accomplished by the chain of 
reasoning given in Table 3, whose geometric and algebraic 
essence is the natural pullback of potentials and forms 
and the natural pushforward of trajectory curves; thus 
the individual steps of Table 3 require only elementary 
manipulations that can be largely automated. 



A. Trajectories are Integral curves 

We first specify the simulation trajectory as an inte- 
gral curve on a metric state-space (Table 3A; henceforth 
"(3A)", etc. for concision). The trajectory increments 
(3A) are specified in terms of one-forms multiplied by 
ltd increments^ dt and dW. Note that stochastic incre- 
ments are expressed by a double-struck "d "; we reserve 
the single-struck "d" for exterior derivative P ' ^^ '. 

In (3B) we discern the dynamical one-form of the 
Schrodinger equation previously described. The familiar 
Dirac braket of quantum mechanics is here regarded as a 
real-valued function on a general Kahlerian state-space; 
we follow Tyurin^^ in calling these functions Berezin sym- 
fto/s'^", although we have retained the physics bra-ket 
notation in token of its familiarity. The symplectic gra- 
dient J*V enters in the Schrodinger dynamics in accord 
with the traverse skiing metaphor of Fig.[l]4_-B. 



C. Lindblad reduction 

That the dynamics of (3A-D) acts to concentrate tra- 
jectories onto lower-dimension manifolds has long been 
noted in the fundamental physics literature. Adler and 
BassP^ have reviewed quantum theories that (in effect) 
concentrate trajectories by gravitationally measuring a 
localized Hamiltonian iJiocai- Especially, Hughston has 
postulated a stochastic state-reduction process— that is 
equivalent to the special case of (3A-C) with q = -ffiocai 
and Sq — {h'^c^/G^y/^/8, where G^ is the gravitational 
constant. Our simulation framework exploits this natural 
concentrative dynamics for practical purposes. 

Thus, ideas from fundamental physics furnish the mid- 
dle panel of a compatible triptych of motivations for 
studying the concentrative one- forms (3A-D): geometers 
conceive them as natural; physicists postulate them as 
fundamental; engineers apply them as agents of concen- 
tration. 



D. The roles of measurement 

Already we can discern the general form of concen- 
tration-and-puUback simulations, and we can inspect the 
dynamic nuclear polarization (DNP) simulation of Fig.p^ 
with good understanding. The output of the simulation 
is a stream of measurements — exactly as in real-world 
experiments. Questions about measurement statistics also 
are well posed, for example: "What is the mean value 
of g"(i)?" These capabilities allow us to duplicate, by 
averaging trajectories and measurements, predictions that 
in Hilbert-space simulation frameworks are achieved by 
density matrix methods. 
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TABLE III. Simulating metric and symplectic quantum dynamics on Kahlerian state-spaces. 

K The multilinearity of tensor network states 
induces block structure in the metric matrix: 



the quantum trajectory equations in differential form 

A Specify trajectory increments via one- 
forms on a state-space manifold JC: 

z{t):t^K (trajectory 
^ ' curve z(t)j 

4,^U-::K^R} (Euclidean 
^ ^^ ■' coordinates) 

(j>''{z{to + St)) = <?i'"(2(tQ)) + / <^"(z{t))dt 
Jto 

d 



9 



,abd(l>\z{t))^Yl 






i-f 



d 



dW 

Ito 



B Drift is the symplectic gradient J*V of a 
Hamiltonian potential (Berezin symbol): 

1-f ^'^ cs 1-f \ Berezin symbol / 1-f 

C Specify Lindblad-lto measurement one-forms 
in terms of Berezin symbols &. their gradients: 






(0I<^) 



measurement 
and noise con- 
centrate states 
by drift toward 
low g-variance 



((l)\qU) 1 continuous 

^^''(^) - ; ' V + ^/"(^) measurement 

('^1'^) \ data g^(£) - 

E\q^{t)q%t')] ^^^SaS{t~t') fxpectation + 

'- ■' ^ -i '■ ' ) Langevin noise 

rto+dt ^ 

■'"' d,„„.,„„ 

_ (^l*^) f {4>\q\4') \ quantum 

p"(8S,)l/2 p^(^|0) jj trajectori 



measurement 
loise drives 
diffusion of 



D The control Hamiltonian H'^(t) is causally 
conditioned upon past measurements q''{t): 

H-{t) = f_^K-{t,t',q-{t'))dt' 

end of trajectory equations, begin state-space pullbacl^ 

E Fullback the Lindblad-lto one-forms onto 
a coordinate manifold !K~G whose image in 
J{ is an immersed simulation state-space X: 

isomorphism quotient map immersion 

3? < — > e — > X — > ^ 

(bijection) (surjection) (injection) 



F Specify real coordinates on X: 

{r":K^K} (coordinates) 

{z"{t) = r''(z(t))} (trajectory) 

{ea = d/dr"} (vector basis) 

{dr^ = V(r")} (dual basis) 
1-f 1-f 
d d 



9at = G^^.^^j (metric G) 

G Specify complex coordinates ^ and holo- 
morphic potentials ip{t;) on !K~C: 

^" if r" + i r" + '1"" ^/^ ; a e 1 , i dim K 

L(t)6K ^''^'/-fi,^'' + <*'-"'^/2; /3el,|dim:K 



ca ig ca^ :p g£f 1 13 (barred index 



{V'''(0:e^C} 



convention) 

(holomorphic 

!/)-potentials) 



H Fullback the Kahler structure of M to C: 
1 d^ 



9alB = g^p 



2 d^^d^!^ 



(^(Olv^(O) 



3a/j = g^i, = 0, 



s 



r , T 1 Mo( 

J inve 



lore- 

trose 

inverse 



I The Kahlerian Lindblad-lto-Stratonovich 
pullback (XLISP) increment is covariant: 

d 



gaf,di^{z{t))^Yl 



fi- 



de 



Strat 



1-f V"? / Strat , 

JU(t)ee 

end of pullback, now factor metric for BLAS computation 

J Introduce tensor network coordinates {■'^*,} 
and specify multilinear ?/;-potentials: 

^(0 = E/C = E/e'®^?'' 



. ® 'C 



algebraic order n- 



\^1 




r 'e 1 




re?i 






pen 




= 




® 




® . 


.® 




^n* 








lc2 






['e;;.,] 



-I- 



r ^ej 1 




r'-e?i 






r '■cr 1 




® 




® . 


.0 




rcl 




["eJ 






r«™J 



te/i = 9aP " 






-■G 



'G 



a/3 



•'■^' G 



{fci}.^^''''} 



2 a^ef ' a^'ef; 



L The product structure of the metric enables 
fast matrix-vector products via BLAS calls: 



[iJG\ = ["c\ (g) ["cj 

metric blocl< vector outer product 



[product 
structure 
of metric 
blocks 



block matrix diagonal matrix 



-^J {*;},{*'('} — °kk' 



\{kl]y '^\{k'l'] 

Comments: (A) Trajectories on a state-space manifold 
Di are stochastic processes whose drifts and diffusions 
are specified by one-forms {^"''., cr"''}. (B) The dynamical 
drift one-form Is the symplectic gradient J*V of a Hamil- 
tonian potential (Berezin symbol function) ((p\H\(p), 
where V is a gradient operator and J* is a complex struc- 
ture on one-forms. (C) The Lindblad-ltO one-forms that 
describe measurement and noise processes are given by 
Berezin symbol functions and their gradients; these act 
to concentrate trajectories onto low-dimension state- 
spaces. (D) The control Hamiltonian is conditioned 
upon on past measurements. 

The Lindblad-lto one-forms given in (C) permit effi- 
cient trajectory integration by dynamically inducing the 
dimension-reducing state-space map (E). The map is 
constructed by introducing real (F) and complex (G) 
coordinates and pulling back holomorphic ^''-potentials 
from J{ to e. This induces a pullback (H) of the Kahler 
structure {J*,G} of D< to C; conversely the image of 
C in Jf is an immersed Kahlerian state-space %. The 
Kahlerian Lindblad-lto-Stratonovich pullback (XLISP) 
increments (I) are covariant and intrinsic to £ and X. 
Concrete trajectories are computed by (numeric) inte- 
gration of DCLISP increments. 

The i/j-potentials are given as tensor network polyno- 
mials (J) whose algebraic order n and Schmidt rank r 
determine the geometry of the state-space X. In turn, 
the induced metric pad on X has a block (K) and prod- 
uct (L) structure such that the matrix-vector product 
gd£, in the !X!LISP increment (I) can be efficiently com- 
puted via basic linear algebra subprograms (BLAS). 



Physically speaking, in the DNP simulations of Fig. [5j 
environmental processes continuously measure and control 
the spin direction, inducing relaxation in accord with spec- 
ified Bloch parameters. High-order quantum correlations 
that otherwise would be costly to simulate are quenched 
by the dynamical flow of trajectories onto low-dimension 
submanifolds. 

A variety of concentration theorems follow from the 
dynamical one-forms of (3A-D)^. In particular, trajec- 
tories can be thermalized and concentrated onto a lower- 
dimension product-space of coherent states by measur- 



ing and controlling all three spin axes simultaneously, 
such that the pulled-back Fokker-Planck equations yield 
positive P -representations of thermal density matrices^. 
Physically, these theorems specify the dynamical one- 
forms of concentration and thermalization. 



E. Quantum trajectory pullback 

Our next task is to pullback the dynamical one-forms 
of (3A-D) onto a reduced-dimension state-space % that is 



A electron-nucleon energy levels B DNP on a rank-1 (classical) state-space C DNP on a rank-2 (Hilbert) state-space 

Tln=T2n = 10 0.6 f ' \ ' ' ] 0.6 
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FIG. 5. Fullback simulation of dynamic nuclear polarization (DNP). (A) The relaxation time constants and equilibrium polarization 
values of an electron-nucleon system are given in terms of Bloch parameters {Tic, r2c, pc} and {Tin, T2n, Pn} (numerical values 
are given in the figure). Switching on the dynamic polarization process at time i = 4 alters the (analytic) mean polarization 
from po = 1/2, Pn = to pc = 3/7, pn = 2/7 (analytic values are drawn as straight lines). (B) Trajectories integrated on a 
rank-1 tensor network state-space have the same dynamical degrees of freedom as a classical two-spin Bloch state-space. Solid 
lines show the fluctuating polarization of the simulated quantum state trajectory; dashed lines show the thermal reservoir's 
(equivalent) continuous measurement of the polarization. Measurement data are low-pass filtered at 6 dB per octave with filter 
poles at 1/(125 Tin) with error bars drawn at ±1 standard deviation of filtered continuous polarization estimate. The simulated 
mean dynamic nuclear polarization is pn — 0.18, which is ~63% of the Hilbert-space value. (C) Trajectories integrated on 
a rank-2 tensor network state-space have the same dynamical degrees of freedom as quantum spins on a Hilbert space. The 
polarization estimates accord with the analytic Hilbert-space values within ±0.02. In trajectory simulations as in experiments, 
lengthy integration time s (7000 Tin in the present case) are required to achieve reasonable statistical accuracy. 



holomorphically immersed in the Hilbert space J{, so that 
"X inherits a Kahler structure by pullback. To accomplish 
this, in (3E-H) we introduce a real coordinate manifold 
3? that is isomorphic to a complex coordinate manifold 
C and we equip C with a set of holomorphic functions 
{ip}- We identify the ■0-functions with (canonical) Hilbert 
coordinates on IK; this induces a natural immersion in 
J{ of the quotient manifold 3C = C/jf/'}; thus a set of 
i/'-functions concretely specifies a Kahler geometry for 
the state-space %. Physically we choose geometries that 
match the concentrative dynamics of the systems we are 
simulating. 

In (31) the Hilbert one-forms of (3A) are specified in 
Stratonovich form; these one- forms pullback covariantljU^ 
and so are geometrically natural. The supporting material 
provides a compendium of Kahlerian-Lindblad-Ito-Stra- 
tonovich pullback (3CLISP) one- forms that is Lindblad- 
complete, that is, sufficient to simulate the trajectories 
of any Lindbladian process (including the Bloch-type 
dynamic polarization process that is simulated in Fig. [sl . 

F. Tensor-network state-spaces 

Now we are ready to choose a class of geometries for 
3C, upon which we will compute the simulations of Fig. [5] 
We choose the tensor network manifolds of (3J). These 
ubiquitouP^ state-spaces are rulcd^'' hyperbolicpS rnani- 
folds; physicists call them matrix product states (MPS). 
The states (3J) are a particular class of MPS whose ma- 
trices are diagonal^, so that r is their Schmidt rank. 
Because multiplication of diagonal matrices is Abelian, 
these states are natural for simulating non-ordered sys- 



tems. Geometrically the states (3J) are a join space^ oi 
multilinear algebraic varieties whose 2"~^ bipartitions all 
have Schmidt rank <r. 

This versatile class of manifolds readily accommodates 
even (non-ordered) classical dynamics, because the pull- 
back of the Hilbert space symplectic form uj onto the 
Schmidt rank-1 product of normalized coherent states 

®fc=r" \Jk,Xk} is simply 



"^Bloch - 



fc=l 



^ jk £abc xl [dx\ ® dxl] 



where jfcG {5,1,!, 



^k'^k^^k 



.} is the fc'th spin's quantum num- 
} is a unit vector on its Bloch 



ber and Xk = { 

sphere; this is the symplectic form of classical spins^and 

(as jfe— )■ 00) of classical test-masseP^. 

A key state-space parameter is the Schmidt rank r; 
as r increases (resp. decreases), the flow becomes 
more quantum (resp. more classical). Example: 
the r — 1 state-space of Fig. [5)3 has dimensionality 
dimjf 3C = (metric) rank [g] = 6; physically these six di- 
mensions are the two Bloch angles of each of n = 2 spins 
plus the two dimensions of a (dynamically irrelevant) com- 
plex amplitude; the resultant dynamics is exactly classical. 
Upon increasing the Schmidt rank to r = 2 in Fig. |5p we 
have dimR 3C = 8 = 2 x 2" = dimR J{; the two added di- 
mensions encode entanglement and the dynamics is fully 
quantum. Thus tensor networks can be tuned to carry 
any required level of quantum entanglement by adapting 
the Schmidt rank. 
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G. Factored metric representations 

To illustrate the computational challenges of state- 
spaces having dimension intermediate between classical 
and Hilbert, we consider a tensor network having or- 
der (say) n — 500 spins and Schmidt rank (say) r = 100. 
This benchmark is selected with a view toward simu- 
lations in imaging and spectroscopy: there are n = 500 
proton spins in 8 nm"^ of water, which is comparable 
to present-day imaging voxels^, and the Schmidt rank 
r = 100 encompasses 2r = 200 degrees of freedom per spin; 
thus the state-space has 100 x the dimensions of classical 
Bloch dynamics, but only an infinitesimal portion of the 
dimffi J{ = 4"=2i™° Hilbert dimensions. 

Computing the trajectory increment d^ amounts to 
solving the JCLISP equations (31) by computing the isomor- 
phism (J)q (called by physicists the raising and lowering 
indices, variously called by mathematicians the canoni- 
cal, metric, natural, or musical isomorphism) that maps 
one-forms onto tangent vectors. A key task in iterative al- 
gorithms for solving matrix equationai^ is computing the 
matrix- vector product (?d^, where g is a 2nr x 2nr Hermi- 
tian matrix; thus naively evaluating g d^ for n — 500 and 
r — 100 requires Ar'^n'^ — lO"'^" multiplications and storage 
elements. 

Seeking efficiency in evaluating ^d^, we find that ten- 
sor network structure induces a metric block structure 
(3K) such that evaluating ^d^ via (3L) requires only 
8r^n = 4-10^ multiplications and 6r^n = 3-10^ storage el- 
ements. The resultant n/2 acceleration and 2n/3 memory 
compression allows the gd^ product for n ~ 500, r = 100 
to be feasibly computed even on small computera^. 

As with Navier-Stokes simulations, the availability of 
0(n) algorithms for matrix- vector products greatly ex- 
tends the range of practical simulation applications. 



H. Physical significance of Schmidt rank 

We have seen that noise levels and Schmidt rank play 
roles in our simulation framework that are broadly analo- 
gous to artificial viscosity and finite-element size in Navier- 
Stokes simulations, that is, they control the trade-off of 
computation cost against simulation accuracy. A strong 
concentration conjecture is that spatially extended, dipole- 
coupled, high-temperature spin systems — such as occur in 
imaging and spectroscopy — are dynamically concentrated 
onto state-spaces of bounded Schmidt rank r; this conjec- 
ture implies that simulating these spin systems is feasible 
with classical resources. Numerical investigations extend- 
ing to 71 = 18 quantum spins on tensor network spaces of 
Schmidt rank up to r = 50 provide some evidence for this 
conjecture^^. 

If further investigations continue to reveal stronger con- 
centration theorems and better numerical algorithms — 
such that puUback remediates (in part) the curse of 
dimensionality — then eventually quantum simulation ca- 
pabilities may parallel those of Navier-Stokes simulation. 



"The linear structure of quantum mechanics is, primarily, 
only a technical convenience" {Ashtekar & Schilling) 
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', "We believe that quantum mechanics [ 

', is in all likelihood exactly linear" {Abrams & Lloyd) ] 

FIG. 6. Hilbert space from a simulationist perspective. Is it 
feasible to follow Ashtekar and Schilling in regarding Hilbert 
space as (|28l p. 61) "primarily, only a technical convenience"? 



V. FIVE NATURAL SIMULATION CHALLENGES 

Now we will suggest five challenges that arise in the 
context of natural simulation. This enumeration is not 
intended to be exhaustive (how could it be?), but rather 
is intended to stimulate readers to envision for themselves 
new challenges and opportunities. 

In drafting these challenges, we were strongly influ- 
enced by DiVincenzo's discussion of dogma and heresy 
in quantum computing^ in particular the DiVincenzo's 
assertion that "We should treat [the quantum] canon 
the way physicists usually treat a canon, by criticizing, 
reinterpreting and creatively flouting it." 

Moreover, a strong engineering tradition is not to criti- 
cize, reinterpret, or creatively flout canons, so much as to 
repurpose them; thus a main focus of our flve challenges 
is to repurpose ideas whose foundations in mathematics 
and science are reasonably well-established, in service of 
practical engineering objectives. 



A. Build quantum-limited microscopes 

In two previous articles (one long22l and one shortP) 
we have reviewed the historical context and the scientific 
rationale of the following challenge: 



Challenge 1 (Quantum-limited spin microscopy). 

Provide first the mathematical and theoretical foundations, 
then the engineering design, then the committed enterprise, 
and finally the practical reality, of spin microscopes capable 
of observing all supramolecular structures, in three dimen- 
sions, with a voxel resolution of one nanometer or less. 



The new capability that natural simulation methods 
bring to this challenge is envisioned to be (with further de- 
velopment) the ability to reliably simulate the entire state- 
space of practical quantum spin microscopes. It is hoped 
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that this capabihty will contribute to an emerging enter- 
prise that is known by many different names — synthetic 
biology, systems biology, and regenerative medicine, to 
name three — that is in reality the single large enterprise 
of achieving a working understanding of our own human 
biology, and our planet's. Foreseeably this will be the 
largest single scientific enterprise of the 21st century; cer- 
tainly it promises to be larger than any previous scientific 
enterprise of any previous century. 



B. Compute fermionic simulations efFiciently 

The natural simulation framework developed in this 
article does not extend to the strongly-coupled fermionic 
systems that characterize much of materials science and 
quantum chemistry. This is not an innate defect, but 
rather reflects mainly the authors' limited assimilation 
(to date) of the vast literature on this subject, as well 
as the practical difficulty of translating methods that 
often are empirical into the language of geometry and 
symplectic dynamics. The following challenge is therefore 
suggested: 



the state-space of quantum mechanics is a vector space. 
Despite efforts, we have to date not found such a natural 
construction for Lindbladian dynamical potentials. We 
therefore pose this challenge: 



Challenge 2 (Efficient fermionic simulations). 

Find efficient, general algorithms for calculating the musical 
isomorphism (that is, T*M-H-TM) on Grassmannian joins 
(in chemists' language, on joins of Slater determinants). 



Challenge [T] and [2] are mutually reinforcing, in the sense 
that advances with respect to either substantially advance 
the other; this is especially evident in the rapidly emerging 
field of synthetic biologji^HUl. 

Challenge |2] — and to a considerable extent all the chal- 
lenges of this section — also is inspired largely by the sem- 
inal work of Ashtekar and Schilling^*, in particular their 
explicit recognition of the desirability of a spin-statistics 
postulate that "refers only to the essential geometric struc- 
tures" and more broadly, their recognition of the need 
for "viable, non-trivial generalizations of quantum kine- 
matics for which even the measurement theory could be 
developed in detail." 



C. Establish Lindbladian naturality 

A conspicuous gap on the overall naturality of our 
simulation framework is that the functional form of the 
Lindbladian potentials of Table |IV| has at present no 
geometrically natural derivation. Instead we derived 
these potentials by a tedious algebraic reduction from 
the Hilbert space expressions presented in Nielsen and 
ChuangiS, guided by the physical ideas of our own recent 
Practical Recipes article^S. More naturally, these Lindbla- 
dian potentials would be constructed solely from intrinsic 
considerations of geometric, dynamical, and informatic 
naturality, that is, without reference to an embedding 
Hilbert space and thus without assuming a priori that 



Challenge 3 (Establish Lindbladian naturality). 

Construct a complete set of Lindbladian stochastic potentials 
solely from, intrinsic considerations of geometric, dynami- 
cal, and informatic naturality (that is, without reference to 
an embedding Hilbert space) and prove that the associated 
communication channels are causal. 



Challenge [3] is intended as a 21st century quantum 
analog to the 19th century challenge of evolving from 
a static Euclidean conception of geometry to a dynami- 
cal Riemannian conception. The latter evolution can be 
broadly conceived as proceeding in four stages: Bowditch's 
1807 discussion of the principles of differential geometry 
is starthng modern (^51 see p. 100); Gauss's celebrated 
Theorema Egregium of 1827 translated these ideas into 
intrinsic terms^ by 1888 Riemann's multidimensional 
generalization of these ideas work was widely known and 
praised^ and yet only in the 20th century was it ex- 
perimentally demonstrated that that the geometry of 
space-time was Riemannian, relativistic, and dynamic. 



D. Demonstrate quantum oracles 

History teaches us that the evolution from Euclidean 
conceptions of geometry to Riemannian conceptions 
aroused passionate debate, and the same is already true of 
quantum geometry. As illustrated in Fig. [61 presently it is 
an open question whether tensor network state-spaces are 
best regarded as a dimensional reduction of Hilbert space 
induced by Lindbladian flow, or conversely, Hilbert space 
might perhaps be merely a dimensional augmentation 
of the Kahlerian state-space of nature, imposed by 20th 
century physicists as a tool for simplifying calculations. 

To focus this debate productively, we therefore pose 
the following challenge: 



Challenge 4 (Demonstrate quantum oracles). 

Experimentally demonstrate any oracle device that cannot 
be simulated in PTIME (that is, with polynomial resources), 
or alternatively, provide experimental evidence that the geo- 
metric dynamics of nature 's quantum state-space obstructs 
such a demonstration. 



Here the idea is to broaden the scope and mitigate 
the technical challenges of quantum information science, 
by embracing some of the ideas of naturality that are 
associated to modern complexity theory. 

To illustrate some of the considerations that arise in 
connection with Challenge |4j let us consider recent work 
benchmarking control methods in 12-qubit spin systemg^. 
Direct simulation of a 12-qubit experiment requires a 
Hilbert space of dimension dimJ{ = 2^^=4096, which 
approaches the practical limit of simulation by density 
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matrix methoda^El. On the other hand, we know from 
the work of Menicucci and CaveaSS that experiments in 
this class in principle can be explained with no quantum 
entanglement whatsoever. The following question thus 
remains unanswered at present: what Schmidt rank is 
required for an accurate trajectory-based simulation of 
a 12-qubit system on a tensor network state-space? The 
simulation framework described in this article is well- 
suited to answering this kind of question. 

We can extend this question to ask, is there any experi- 
ment we can accomplish, whose accurate simulation prov- 
ably is infeasible with PTIME computational resources? 

In our view, such an experiment would be comparable in 
scientific signifiance to Galileo's (possibly legendary) ball- 
drop from the Tower of Pisa. Challenge |4] is framed with 
a view toward helping to catalyze such an experiment. 

It is tempting to frame these issues as a 21st century 
struggle between Boeotian believers in a Hilbert state- 
space^, versus skeptics of quantum computing^. A more 
traditional and sober view, though, is that a sustained 
tension between quantum Boeotians and quantum skeptics 
is desirable as a stimulus to progress in quantum infor- 
mation science. As ever, strong experiments analyzed in 
the light of well-posed theories have the final say. 



E. Create natural field theories 

The following question originated in discussions within 
our Quantum Systems Engineering (QSE) Group: Why 
is our universe friendly to simulation? 



Challenge 5 (Naturality in field theory). 

Construct a relativistic field theory whose near-field dynam- 
ics exhibits Lindbladian naturality, whose communication 
channels are provably causal, and whose quantum and space- 
time structures (that is, whose Kahlenan and Riemannian 
structures) are dynamically coupled. 



As with the preceding challenges, Challenge [5] too has 
a dual aspect. Even if it should happen that the quan- 
tum state-space of Nature is non-dynamical and perfectly 
Hilbert, there still may be substantial computational ad- 
vantages attendant to regarding quantum state-space as 
dynamical and Kahlerian; here the point is that efficient 
techniques for simulating quantum field dynamics has 
many practical uses in optics, radar, and condensed mat- 
ter physics. 

On the other hand. Challenge [5] encourages us to con- 
ceive of field theory along lines that differ substantially 
from 20th century conceptions, in which information- 
theoretic and complexity-theoretic considerations are sim- 
ilarly important as dynamical considerations. We can 
envision, for example, an outcome along the lines "M- 
theory is the unique, causally separable, relativistically 
invariant quantum field theory that can be simulated with 
PTIME computational resources." 



VI. CONCLUSIONS 

The elements of naturality that we have described have 
allowed us to naturally specify a dynamical simulation 
framework that encompasses Hamiltonian, Lindbladian, 
and thermostatic flows on classical and quantum state- 
spaces. 



A. Educational considerations 

Naturality is purchased at significant price: our de- 
scription employs approximately one hundred abstract 
mathematical terms over-and-beyond those terms that 
typically are taught in the context of an undergradu- 
ate education in mathematics, engineering, and science. 
Excellent textbooks exist for learning these elements of 
naturality — such the dozen canonical texts^Hl^ suggested 
in Sec. |IE| — but at present there is no textbook that 
presents these naturality concepts in the unified context 
of practical simulation; perhaps someday textbooks on 
simulation may be written taking this perspective. 

Once the canons of naturality are mastered, the sim- 
ulation framework presented in this article itself flows 
naturally from elementary considerations of geometry, 
dynamics, and information theory. Only a handful of 
equations need be grasped, and key proofs are associated 
to to projective diagrams whose form is easily recalled 
(see Theorems 1-3 of Sec. Im] and Tables IM!]) . 



B. Unleashing capabilities from tethers 

The flve natural simulation challenges whose technical 
aspects are discussed in the preceding section (Section [y]) 
can be understood as being collectively focussed on repur- 
posing existing capabilities and creating new capabilities 
that unleash us from present-day tethers, as follows. 

Challenge [l} Quantum-limited spin microscopy 

The spin microscopy challenge seeks to unleash us from 
the tether of pursuing structural and synthetic biology, 
and concomitantly regenerative biomedicine, primarily by 
experimental methods. As our recent overview of spin 
microscopy concludes ii' 

We are tantalized by a vision of medical prac- 
tice becoming fully curative and regenerative. 
We are frustrated — as the generation of von 
Neumann and Wiener was frustrated — by the 
limitations of our present tools. We desire — as 
Feynman famously desirecP^ — to "just look at 
the thing" . And we plan — as every previous 
generation has planned — for these aspirations 
to become realities. 

The present article is thus the third and flnal arti- 
cle in a tryptycfpESl ii^sX focusses on von Neumann's 
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challenge^ to consider whether atomic-resolution biolog- 
ical microscopy might be achieved "by developments of 
which we can already foresee the character, the caliber, 
and the duration. And are the latter two not excessive 
and impractical?" 

The natural simulation framework presented in this 
article was conceived, first and foremost, to help 21st 
century mathematicians, scientists, and engineers to meet 
von Neumann's biomicroscopy challenge, and thereby to 
dynamically augment of knowledge of healing, by creating 
a balanced surge of new biomedical research capabilities. 

Challenge [2| Efficient Fermionic simulation 

The Fermionic simulation challenge seeks to unleash 
condensed matter physics from the tether of Fermi liquid 
models, by providing efficient numerical tools for exploring 
condensed matter physics beyond Fermi-type approxima- 
tions. This challenge seeks to bring to materials science a 
predictive power comparable in scope to the observational 
power that quantum spin microscopy seeks to bring to 
biomedical science. From a practical point of view, this 
augmentation of materials science capability will play a 
strategically vital role in the 21st century, as spurred by 
urgent necessity, humanity seeks new ways to free itself 
from age-old tethers like carbon-burning and resource 
shortfalls. 

Challenge [3} Establish Lindbladian naturality 

The Lindbladian naturality challenge seeks to expand 
the predictive power of quantum physics by unleashing 
it from the tether of Hilbert space. No matter whether 
the state-space of Nature is a Hilbert space or not, by 
recognizing that Lindbladian flow is generically dimension- 
reducing, we gain immense new computational powers 
and design insights in the quantum realm. 



Challenge |4} Demonstrate quantum oracles 

The quantum oracle challenge seeks to expand the realm 
of quantum information science by unleashing quantum 
information processing from the tether of quantum com- 
puting. The key recognition is that the quantum oracle 
challenge is mathematically broader and potentially easier 
of experimental achievement than quantum computing, 
and yet is similarly seminal in its influence on our under- 
standing of the computational aspects of our universe. 

Challenge [Sj Naturality in field theory 

This mathematically toughest of our five challenges 
seeks to expand the capabilities of field theory by unleash- 
ing it from the tether of dynamical dominance. The key 
recognition is that informatic flow is equally central to 
dynamical flow, and that the quantum and space-time 
structures of nature's state-space are (potentially) natu- 
rally dynamic and mutually coupled. 

We close by emphasizing that the above five challenges 
(and many more that might be conceived) are closely 
dovetailed, in the sense that progress on any one of them 
acts to speed the pace and retire the risks of progress 
toward meeting all of them. 
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Appendix A: Lindbladian recipes 

Given the general Block parameters of a quantum spin 
system, this supplement provides explicit expressions for 
that system's Kahlerian Lindblad-Ito-Stratonovich pull- 
back (DCLISP) one-forms (31). These expressions suffice 
to code the simulations of the main article, and more 
broadly, they suffice (in principle) to simulate all quan- 
tum systems whose dynamics are Lindbladian, including 
classical systems that respect the standard quantum limits 
(see Sees. 3.2.8 and 3.4 of^. 

Three key elements of geometric simulation methods 
are summarized in this appendix: first, a Bloch-complete 
tabulation of Lindblad-Ito generators (4A-E), second, a 
general covariant (Stratonovich) form of the 3CLISP incre- 
ment (4G-I), and third, a discussion of practical implica- 
tions of Kahlerian metric factoring. 



1. Completely positive linear maps 

To begin, it is well known (fHl see Sec. 8.4.1) that 
the most general infinitesimal generator of a completely 
positive linear map on a density matrix p (physically, a 
map that preserves the positivity of probabilities) is of 
the following Lindbladian form: 



Lip) = -t[H,p]/h 

n 

+ E 



fe=i 



2Sa, 



LkpLk ~ hP^kLk - \L,,L 



kP 



(Al) 



where H \s a, (Hermitian) Hamiltonian operator, the 
{i^; k € 1,.. .n} are a set of (nonHermitian) Lindbla- 
dian generators, and a superscript " H " indicates Her- 
mitian conjugation. Without loss of generality we take 
Lk ~ Ak + iBk , with Ak and Bk being Hermitian opera- 
tors and we introduce a (real, positive) scalar normalizing 
factor 5'^^ that we will identify with a (one-sided) mea- 
surement noise spectral density having units of A^ ■ time. 
Note that the normalizing factor Sa^. in Eq. [XT allows 
the operators Lk to be assigned any desired physical 
units (for example, position, linear momentum, angular 
momentum) . 



2. Quantum operation unravellings 

We focus our attention upon the Lindbladian map that 
is generated by a single (arbitrary) Lindbladian opera- 
tor L. Following Nielsen and Chuang (|141 chs. 2 and 8) 
we first exhibit a stochastic unravelling of a simulated 
wave function ip{t) that is specified in terms of quantum 
operations M~^ and M~ by 



V'(t + St) 



M^ijj{t)/ ^Jp+ with probability p+, 
M^ip{t)/ ^Jp'' with probability p" 



where the probabilities p are given by 

p± - (M±^(t)|M±^(t)), (A3) 

and the form of the measurement operations M^ is 

M± =exp ( - i{AB + BA)5tl{ASA)) 
xejip{±iB^ySt/{2SA)) 



X [cos {A^y St /{2S a)) 

±sm{A^St/{2SA))]/V2. (A4) 

By the physic al re asoning oi^ (Sec. 3.3 and Fig. 9), we 
recognize Eq. (A4) as a synoptic unravelling, with A as 



the measured operator, B as the feedback Hamiltonian 
operator, and the term oc {AB + BA) as a Stark shift 
Hamiltonian. 

The time-evolution of p then is given in terms of the 
M± by 

p{t + St) - p{t) = M+p{t){M+)" 

+ M-pit){M-)^~p{t). (A5) 



Substituting Eq. ( [A2p3| in Eq. ( [A5| and expanding 
in VSt immediately establishes the equivalence of the 
quantum operation unravelling to the Lindbladian map: 



L(p) = lim 

1 
~2S^ 



p{t + St) - p{t) 
St 

[Lp{t)L'' 



y{t)L^L 



(A6) 
\L»Lp{t)\. (A7) 



(A2) 



It is important to appreciate that the above Lindbla- 
dian map does not uniquely define the generator L\ it 
is invariant under L — > iL, which corresponds to A^ B 
and B — )■ —A. The operator A is privileged in the unrav- 
elling by associating it with a measurement process; in 
consequence it is solely A that will appear in the crucial 
(for purposes of concentration) variance-reducing drift 
one-forms of (4G) and (41). 



3. Finite-polarization Bloch equations 

Recipes for elaborating the above Lindbladian map 
into general-purpose recipes for simulating spin systems 
described by Bloch-type relaxation are well covered in the 
literature (Sees. 3.4 & 4. 1(23), and (4B-E) summarize the 
needed relations. 

In brief, Bloch-type relaxation to thermal equilibrium 
is unravelled as a process of measurement and control, 
in which the environment monitors the spin direction 
and applies control feedback necessary to establish the 
polarization is attendant to thermal equilibrium. 

More specifically, the operation representation of (4B) 
unravels the thermalizing process as the sum of three 
measurement-and-control operations, in which the ex- 
pectation value of each of spin operators {sx,Sy,Sz\ is 
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TABLE IV. Recipes for the simulation of Lindbladian quantum dynamics on Kahlerian state-spaces. 



Specify the Lindblad & Bloch dynamics 

A Specify a general Lindblad superoperator £ 
in terms of generators {Lk = A^ + iBk}: 

C{p) ^ dp{t)/dt ^ -i[H,p]/h 

+ E ^ [L^pL'i - \PL1L, - \niL,p] 
k=i ^^■■^'' 

B Specify a thermal Lindblad superoperator £'*■ 
whose equilibrium state is p"' oc (!xp(— /Ss^) 

£), = [L.pn', - \pL",L, - \L",L,p]/{2S^) 
+ [LyPLi - kpL"yLy - \L"^Lyp]l{2Sy) 
+ [L^pL^^ - \pL"^L, - \L"^L,p]/{2S,) 

The thermal Lindblad generators are given 
in terms of spin-j operators {sj,,,Sj,,,Sj} by: 

Lx = *i ^ iaSy ;i-axis measure-and-control 

Ly = Sy + ias.j. j/-axis measure-and-control 

i^ = s^ 2-axis measure only 

a = tanh(/3/4) thermal control gain 

C Specify Bloch parameters {Ti,T2,po}: 



d_ 
dt 



tr s^p{t) 
trsyp(t) 
tr Szp{t) 



tr Sxp(t)/T2 

trsyp(t)/T2 

tr(s,-jpo)pit)/Ti 



end Lindblad & Bloch dynamics: specify conversions 

D Specify the Bloch parameters {po.Ti,r2} in 
terms of Lindblad parameters {a, Sx = 5^, 5^}: 

"=((l-Po)'/'-l)/Po 

s, = (4/r2 - 2/ri)-i 



E Specify Lindblad parameters {a,S';r = S^, 5^} 
in terms of Bloch parameters {po:Ti,T2}: 

po = -2/(a + 1/a) 

Ti = 2SJ{l + a^) 

T2 = 45,5,/(5,, + 5,(l + Q^)) 

end conversions, specify the Lindblad-lto dynamical forms 

F Unravel L (see A) in terms of concentrative 
Lindblad-lto one-forms {o", /(} as specified 
by Berezin functions and their gradients: 

L = A + iB (A, B both Hermitian) 

/b(-) : 3C ^ K i=' (01 • 10) 
then the diffusion one-form a is 

7b(A) 



Berezin's 

symbol 
functional 



{85a)i/^ 






MI) 



where / is the identity operator. 

G Measurement of operator A with PSD Sa 
induces Hamiltonian noise with PSD Sh'- 

Sh = H^A^Sa 

H Specify the Lindblad-lto drift one-form p 
(but p pullback C <- 3C <- J{ not covariant): 



/-* = 



85. 



-V(/b(A^: 
'MA) 



MA)-" 

MI) 

,/B(/)/*^^^(^) 
J*\7MAB-BA) 

■r^MH) 



-^MB') 



I Specify the Stratonovich drift one-form /( 
(now p. pullback 6 <— DC <— 3i js covariant): 



1 

''&Sa 



-2/b(/)V 
\,M^ 

MI) 



MA^) 
MI) 

r^MB) 



MA)' 
MI)' 



t (T) r.^( MAB-BA) \ 



-.rVfB{AB+BA) 



2h 



2h 



Comments: (A) The Lindblad superoperator L 
acts linearly upon a density matrix p; here "H" 
is Hermitian conjugation. (B) Bloch-type relax- 
ation is modeled in terms of three Lindblad gen- 
erators {Lx,Ly,Lz} given in terms of spin oper- 
ators [sx,Sy]—iSs and a control gain a. Bloch 
relaxation is induced by three-axis measurements 
whose (one-sided) power spectra! densities (PSDs) 
are {Sx, Sy, Sz}- (C) Relaxation of the density ma- 
trix in terms of Bloch parameters. (D) The Bloch 
parameters in terms of the Lindblad parameters. 
(E) The Lindblad parameters in terms of the Bloch 
parameters. (F) The diffusion one-form associated 
with the synoptic unravelling of an arbitrary Lindblad 
generator L. By definition /b(^) ^^ (01-^ 1^) is the 
Berezin function of a Hermitian operator X . Physi- 
cally A is a measured operator and B is a feedback 
operator; thus A induces Hilbert backaction via a 
metric (V) gradient white B induces Hamiltonian 
flow via a symplectic ( J*V) gradient. (G) The invari- 
ance of the superoperator L under L—^iL ensures 
that continuous observation of an operator A is op- 
erationally equivalent to application of a stochastic 
Hamiltonian. (H) The Ito drift one-form for the syn- 
optic unravelling of a Lindblad generator L = A + iB. 
(I) The Stratonovich drift induced by the Ito drift. 



continuously measured by processes having spectral den- 
sities {Sx,Sy,Sz} (one-sided), and feedback control is 
applied with gain a. To assign physical values to the four 
thermalizing parameters {5*2:, S'j,, S'2,a}, we first specify 
Sx = Sy (transverse isotropy) . Then by a purely algebraic 
computation that depends only upon the spin commuta- 
tion relations, we find that setting a = tanh(/3/4) ensures 
£(exp(— /^s^)) = 0, i.e., the thermal density operator 
Pth = exp(— /Js^) is a fixed point of the thermalizing 
Lindbladian map. 

In this picture the control gain — the parameter a in 
(4B-E) — sets the temperature, with optimal gain of a = 1 
corresponding to zero temperature, and lesser gains (or 
larger gains, see^ corresponding to finite temperature. 

This physical picture explains why we cannot simu- 
late quantum computations via geometric methods: the 
covert observation of the thermal environment destroys 
the high-order correlations upon which quantum compu- 
tation depends, and which otherwise would be costly to 
simulate. Present-day designs for quantum computing 
rely on an external supply of ancilla bits to error-correct 
these high-order correlations; these bits represent a steady 
increase in the dimension of the Hilbert space that simi- 



larly is not accommodated in the tensor network manifold 
upon which 3CLISP increments are pulled-back. 

The relations of (4B-E) are valid for arbitrary spin-j, 
and in particular the limit j' — > 00 describes the motion 
of test-masses; this limit too is discussed in-depth irP3 
(Sees. 3.4 & 4.1). 



4. Integrating stochastic equations 

By a straightforward calculation, we can check that the 
synoptic unravelling of Eqs. (A2-A5) is represented by 
the Ito increment one-forms of (4F-H). However, in order 



to pullback these one-forms onto a lower-dimension %- 
manifold and then numerically integrate the trajectory, we 
must calculate a correction term called the Stratonovich 
correction. 

We now derive this correction. We introduce Kahlerian 
indices as in the main article such that the Hilbert metric 
g'K is given in terms of the Kahler potential k(i/)) — 
{i^\i^)/2 = /b(/)/2 by 



[5m] ab = Gji 



d d 



1 



92 



di}}''' dtp'' J 2 dip^'dtp'' 



/b(/), (A8) 
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where /b is Berezin's symbol functional (see (4F) for defi- 
nition) and / is the identity operator on the Hilbert space 
3{. 

A technical point is that our definition of Berezin's 
symbol functional is less restrictive than that of most 
authors in geometric quantum mechanics; for us /b is 
simply a linear map from Hermitian operators to real- 
valued state-space functions. This reflects the fact that in 
our formalism Hamiltonian flow is a symplectomorphism 
(the Lie derivative of the symplectic form is zero because 
the pullback of a closed symplectic form is closed) but 
not necessarily a (metric) isomorphism (a flow that leaves 
the metric form invariant; compare Proposition 1.3 oi^SI). 

We do not require that pulled-back Hamiltonian flow 
be a (metric) isomorphism because we know from prior 
theoretical work^'' that tensor network state-spaces have 
nonuniform Riemannian curvature structure, such that 
Hamiltonian trajectories span regions of large and small 
state-space curvature. Furthermore, as discussed in the 
main article, we know that symplectomorphic flow is a 
strong enough condition to ensure thermodynamically 
reasonable simulation physics. 

As in the main article we define the raised tensor com- 
ponents [ffjt]"^ to be the (unique) matrix pseudo-inverse 
of [g:H]ab- 



[g^r [9^]bc [5M 



\cd 



[5M 



\ad 



(A9) 



The metric gx on the pullback Kahler state-space % is 

(AlO) 



r 1 ^V-'r 1 dtp" 



Now we can define a (local) projection matrix [Pjcl^'h in 
terms of the metric on our state-space: 

[Poc]^ = ^^9oc]-^[9:KU. (Alt) 

We consider the Hilbert-space Ito increment from (4A) of 
the main article: 



dV'" = [5Jf 



]ab 



A^l^^ld* 



-'^^''^^ 



(A12) 



Applying Pjc to Eq. ( A12 1 and using the chain rule, we ob- 



tain a covariant expression for the Ito increment projected 
onto the tangent state-space 



[p^nd^y 



b _ ^^^ r„„.i™« 



ae' 



m 



I'l^^'^t 



'y\g^)dW 



(A13) 



Now we let the £, functions be processes with Ito incre- 
ments {/z^jCr"} that we wish to solve for 



dC = fi'^dt + a°'dW. 



(A14) 



Ito's Lemma then gives the ?/;-increments in terms of 
{/^°,cr"} as 

f)al,b 1 f)2nl.b 

dV'''^^ [^iPdt + aPdW] + -—^a'^a'dt. (A15) 



Applying Px yields 

dtp" 



[Pxndr 



5e 



goc\ 



[9oc 
9^V 



np 
2„;,6 



^lPdt + aPdW 



lM-'^[9.U^^^a^a^dt 



(A16) 



We equate the right-hand-sides of Eqs. ( A13) a nd (A16 l, 
and switch to Kahlerian indices. Recalling Eq. (AlO I, we 



readily verify that the resulting equation is satisfied by 
the following ^-coordinate Ito drift and diffusion: 



." = [^^r-(|^), 



M 



[9x 



a/3 



f d \ ld[gx]^ , 



(A17) 
(A18) 



This expression has the computationally desirable prop- 
erty of being wholly intrinsic to the !K-manifold. However 
we notice that the pullback C <— 3C <— IK has induced a 
noncovariant correction to the Ito drift (because the fac- 
tor d[gx]j3^/d£,^ in Eq. (A18l does not transform as a 
tensor) . 

We see that Ito's Lemma generically induces noncovari- 
ant drifts whenever Ito increments are pulled-back onto 
nonEuclidean manifolds. This implies that Ito drift one- 
forms are not natural geometric objects; in the following 
section this will motivate us to define Stratonovich drift 
one-forms that are geometrically natural. 

However, Ito increments are quite convenient for numer- 
ical computations, and the dynamic nuclear polarization 
(DNP) trajectories reported in this article were in fact 
computed by integrating Ito increments. Remark: the 
high-precision (±0.02 in polarization) agreement between 
density matrix and geometric trajectory calculations that 
is seen in (3C) required the inclusion of the correction 
term in Eq. (A18); this finding was contrary to our early 
expectations (or hopes) that the correction term might 
be dynamically negligible. 

More broadly, the expressions collected in Table 3 of the 
main article, together with the compendium of Table 4, 
plus the above drift correction of Eq. (A18), suffice to 



numerically simulate any spin system whose symplectic 
dynamics are Hamiltonian, whose metric dynamics are 
Lindbladian, and whose relaxation parameters are Bloch- 
type. 



5. Ito-to-Stratonovich conversion 

A well known, formally equivalent, and geometrically 
covariant alternative to Ito integration is Stratonovich 
integration, whose Kahlerian drift and diffusion we now 
derive. 

In essence Stratonovich methods embrace an alternative 
stochastic integration scheme that eliminates the nonco- 
variant correction term in Eqs. (A15I and (A18|, such 



that the drift and diffusion one-forms pullback 6 <— 3C ^ IK 
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covariantly, and thus retain a natural geometric interpre- 
tation. This geometric naturalness is obtained at the 
expense of a more complicated, nonMarkovian definition 
of stochastic integration; this trade-off is the subject of a 
large literature. In general Ito and Stratonovich integra- 
tion methods both are widely used, and the decision of 
which to use is made on a case-by-case basis. 

On the Hilbert manifold !K the Stratonovich drift one- 
form fi is given covariantly in terms of the Ito drift one- 
form /i by a well known conversion relation: 



d 

dip"- 



d 

dip'' 

1 r ,bc f d \ d 



^V". 



(A19) 



Switching to Kahlerian indices, but still working on IK, 
we obtain 



d 



dip" 



d 



dip' 



[9:Kra 



d \ d 



dip^J dipP \ \dip" 



d 



[g^ra 



d 



d 



dip-'J dip^ \ \dip 



d 



(A20) 



Substituting the Ito drift and diffusion from (4F) and (4H), 
carrying through the indicated derivatives, and expressing 
the result in terms of Berezin functions and gradient one- 
forms, yields the Stratonovich drift one-forms of (41). 

We find that the Stratonovich drift of F (41) is variance- 
reducing and thus concentrative, just like the Ito drift of 
(4H). 

The chief practical advantage of Stratonovich drifts 
for simulation purposes is that they pullback covariantly 
onto the reduced-dimension state-space 3C; in particular 
the noncovariant (and computationally costly) correction 
terms of the preceding section are not induced by the 
pullback. 



6. Kahlerian metric factoring 

In our own research we have not as yet numerically ex- 
perimented with Stratonovich integration; however Kloe- 
den and Platen's textbook^^ suggests specific lines of 
inquiry for stochastic systems in this class. For example, 
the state-space flows induced by the diffusive one-forms 
of (4F) commute in the stochastic sense if and only if 
their respective Lindbladian generators commute. Numer- 
ical stochastic integration — both Ito and Stratonovich — 
can be made more efficient whenever all (or most) diffu- 
sions commute; this diffusion structure occurs naturally 
in Bloch-parameter spin models. 



As was discussed in the main article, computing the 
gradient-to- vector isomorphism for a system of n = 500 
spin- 1/2 particles, pulled-back onto a rank r = 100 tensor 
network state-space, requires the solution of 2nr — 10^ 
simultaneous equations. Such systems are too l arge for 
direct numerical solution, and it is well knowrpSEMl 
that iterative methods can lead to O(n)-fold speed-up 
whenever the associated matrix-vector product can be 
evaluated in 0(n) operations. For systems having n of 
order 10^ or larger, this can be the difference between 
feasible and infeasible simulation; for this reason iterative 
methods are almost universally employed in large-scale 
classical simulations. 

Iterative methods, generally speaking, are optimized 
on a case-by-case basis; thus it is reasonable to anticipate 
that the development of iterative methods for integrat- 
ing simulation trajectories will proceed by a blend of 
fundamental theory guided by practical computational 
experience. 



7. Observations upon trajectory concentration 

Let us examine the expressions for 3CLISP diffusion (4F) 
and Stratonovich drift (41) with a view toward identifying 
aspects of trajectory concentration that are understood 
versus not understood, in both mathematical and physical 
terms. 

Turning our attention first to the diffusion one-forms 
(4F), it is apparent that the symplectic gradient one-form 
ex J*V/b(-B) generates the Hamiltonian fiow of the feed- 
back noise, while the metric gradient one- form ex V/b(^) 
generates the Hilbert backaction that is associated with 
the measurement of A. It is this feedback that is dynami- 
cally responsible for familiar quantum phenomena such 
as the Stern-Gerlach effect, for example (see^^ (Sec. 3.2)). 

Turning our attention now to the drift one-forms (41), 
the sole metric gradient one-form is manifestly concen- 
trative; it is responsible for the existence of the P- 
representations described in Sec. 3.2 oP^l, for example. 
The symplectic gradient ex J*y f^{B) supplies the coher- 
ent feedback associated with measurement of A. The com- 
mutator term AB -\- BA (which functions as a Stark shift 
Hamiltonian) was inserted at the beginning (Eq. A4 ) solely 
to obtain the conventional Lindbladian map (Eq. Al). 
The anticommutator term AB — BA presumably has geo- 



metric significance, but this significance is opaque (to us 
at least). 

The unravelling we have been discussing is synoptic, 
and thus the preceding discussion is informed by a physi- 
cal picture in which A is measured and B specifies feed- 
back. However, the synoptic representation is not unique; 
there exists a unitary invariance — called the Theorema 
Dilectum (theorem of choice)p21 — which guarantees that 
a class of informatically equivalent symplectic and met- 
ric one-forms can be unravelled from any given starting 
Lindbladian map. This theorem is manifest in the invari- 



ance of the Lindbladian generator L of Eq. ( Al ) under 
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the mapping induced by L^ i-)- ^ UkjLj, where U is an 
arbitrary unitary matrix. The physical and informatic 
significance of the [/-matrix is discussed irP21 in essence 
L'^-equivalence is necessary to the causal separability of 
quantum measurements. 

Until we can survey the entire [/-equivalence class of 
0<!LISP one-forms with mathematical equanimity, it can- 
not be said that we have a satisfactory understanding of 
the geometric aspects of quantum dynamics. Similarly, 
until we can select — or better, design — particular unravel- 
lings that are well-suited to efficient simulation, it cannot 
be said that we possess even a working grasp of quantum 
simulation's computational complexity. 

This uncertainty extends to simulation's close cousin, 
optimization. Suppose for example that we have a clas- 
sical annealing algorithm for solving (inefficiently) an 
NP-complete problem like 3-SAT. If we view the 3-SAT 
bits as qubits, and the annealing process as a quantum 
annealing process, then we can ask questions like: "Bear- 
ing in mind the gauge-like invariance of the Theorema 
Dilectum that is associated to Lindbladian dynamic^^, 
what quantum noise model yields the most effective an- 
nealing on a Schmidt rank-r state-space?" Section 4.6 
oP^ discusses these questions, but informatically speaking, 
we know very little about optimizing flow on Kahlerian 
state-spaces. 
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